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Abstract 

We investigate the effect of heavy ffavor contributions to vector boson {V = 7, Z, 
W) production which is described by the Drell-Yan mechanism. All reactions with 
bottom and top quarks {Qi = b, t) in the final state, like gi + g2 V + Qi + Q2 
and g + g V + Qi + Q2, are considered. This study also includes the virtual 
contributions containing heavy flavor loops which were not taken into account earlier 
in the literature. Our analysis reveals that the above corrections to the Drell-Yan 
cross section are very small. Only at energies characteristic for the LHC they are of 
the same order of magnitude as the order QCD contributions due to light quark 
and gluon subprocesses calculated earlier in the literature. 



1 Introduction 



The production cross sections for the electroweak vector bosons W and Z in hadron- 
hadron coUisions as described by the Drell-Yan (DY) [1] mechanism provides us with 
a beautiful test of perturbative QCD. One of the reasons is that the total cross section 
can be calculated in next-to-next-to leading order in perturbative QCD [2, 3] a result 
which is very hard to achieve for other processes in hadron-hadron collisions. 
Another advantage of this process is that on the Born level it is purely electroweak in 
origin for which the theory is in an excellent shape. Hence each deviation of the cross 
section from the Born approximation can be attributed to QCD effects. Therefore 
the Drcll-Yan (DY) process belongs to the same class as deep inelastic Icpton hadron 
scattering and electron-positron collisions where QCD corrections can be measured 
with much higher accuracy than in other reactions. 

Besides the total cross section for vector boson production and the cross section 
da /dm where m denotes the lepton pair invariant mass order corrections due to 
soft plus virtual gluon contributions have been calculated in [4] to the differential 
distributions d^a/dmdxp and d'^a/dmdy. Here the quantities xp and y denote the 
fraction of the longitudinal momentum of the lepton pair with respect to the center 
of mass (CM) momentum and the rapidity respectively. 

Furthermore one has computed the full order as correction to the single vector bo- 
son inclusive cross sections d^a /dmdxpdpT or d^a/dmdydpT [7] where pt denotes 
the transverse momentum of the vector boson. All the above calculations have been 
performed under the assumption that the quarks, appearing in the partonic subpro- 
cesses contributing to the DY process, are massless. This is a reasonable assumption 
for the light quarks u, d, and s and even for c since the masses of the W and Z 
are large compared with the masses of the above quarks. However this assumption is 
doubtful for the bottom quark and it is certainly wrong for the top quark since recent 
experiments [9] indicate that rrit > My\r, Mz- Therefore we cannot neglect the masses 
of the bottom and top in the final state of the partonic subprocesses in particular if 
the collider CM energy is small like in the case of the SppS {^/s — 0.63 TeV) or the 
TEVATRON (y^ = 1.8 TeV). Maybe for some partonic subprocesses one can apply 
the zero mass approximation for the bottom quark at LHC (y^ = 16 TeV) as we will 
show later on. Contrary to the final state we will omit the bottom and top quarks in 
the initial state because we assume that the bottom and top densities in the hadron 
are negligibly small. 

The calculation of the contribution of heavy flavors to the DY cross section has been 
performed for Z-production in [10]. It contained all one-loop and two-loop corrections 
which are characteristic of Z-production but do not contribute to 14^-production or 
processes with a photon in the intermediate state. They are characterized by the tri- 
angular heavy flavor-loop insertions containing the Adler-Bell-Jackiw anomaly which 
has to cancel while adding top and bottom loops. This work was extended in [11] by 
including the interference terms originating from the final and initial state radiation 
of the vector boson in the process q + q^Z + Q + Q (with Q = b, t). All other 
production mechanisms, which also show up for ^/-production and processes with a 
photon in the intermediate state, have not been considered yet. The contributions to 
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Z-production considered in [10, 11] all show up on the order level and amount to 
about 0.1 percent of the Born approximation which means that they are experimen- 
tally unobservable. 

In this paper we want to complete the above analysis by including all remaining 
Feynman graphs which also contribute to lepton pair production with a W or a 
photon in the intermediate state. Apart of some additional two-loop graphs they 
contain the contribution of the partonic subproccsscs qi + q2 V + Qi + Q2 and 
g + g ^ V + Qi + Q2 with V = ^, Z , W and Qi = b, t. Like the corrections discussed 
in [10, 11] they all contribute to the DY cross section on the order a1 level. 
This paper will be organized as follows. In section 2 we present the partonic cross 
sections corresponding to the subprocesses which contribute to heavy flavors plus vec- 
tor boson production. Furthermore we show that in the case the vector boson mass 
becomes much larger than the heavy flavor mass one can find explicit analytic results. 
In section 3 we compute the heavy flavor contributions to vector boson production at 
current and future hadron-hadron coUiders. In particular a comparison will be made 
between the order corrections due to light partons (quarks and gluons) , calculated 
previously in [2, 3], and the contributions due to bottom and top quarks presented 
in this paper. Finally in Appendix A we give an explicit formula for the two-loop 
vertex correction which was not computed in [10, 11]. 
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2 The order corrections to the DY cross section 
due to heavy flavor production 

In this section we present the partonic cross sections of heavy flavor production 
contributing to the Drell-Yan (DY) process which is given by 

H1 + H2 y + 

U£,+4, (2-1) 

where Hi, H2 denote the incoming hadrons and V stands for one of the vector bosons 
of the standard electroweak model (7, Z or W) which subsequently decays in the 
lepton pair £1, £2. 

The heavy flavors in the final state are given by Qi and Q2 respectively and the symbol 
"X" denotes any inclusive final hadronic state. In lowest order of the electroweak and 
strong coupling constants the above reaction receives contributions of the following 
parton subprocesses 

i{ki) +j{k2) ^ V{q) + gi(pi) + Q2{P2), (2.2) 

with i, j — q, q, g. Here q stands for the light quarks given by u, d, s, and c whereas 

the heavy qiiarks are represented by t and b. Notice that in this paper we study 
heavy flavor production at large hadron collider energies so that the charm can be 
treated as a light quark which mass can be neglected. In addition to reaction (2.2) 
we include the virtual corrections due to heavy flavor loops which contribute to the 
subprocesses 

tih)+jik2) ^ Viq), (2.3) 

i{h)+j{k2) ^ Viq) + l{ks), (2.4) 

with I — q, q, g- The most of these virtual corrections reveal the presence of the 
triangle fermion loops giving rise to the well known Adler-Bell-Jackiw anomaly which 

has to cancel between the top and bottom contributions. They have been treated 
in [10, 11] in the case of V = Z. Here below we will add the contributions due to 
gluon self energies containing the heavy quark loop which also appear when V is 
represented by 7 and W. 

In the subsequent part of this paper we are interested in the following parton cross 
section 



do-V- CP- 



T av{Q\ Ml) m{f, Q\ ml ml), f = ^, (2.5) 



receiving contributions from the reactions (2.2)-(2.4) {i, j = q, q, g). The quantity 
(Ty represents the pointlike DY cross section and the kinematical variables V§ {s — 
(ki + k2)^) and y/Q^ stand for the CM. energy and the lepton pair invariant mass 
respectively. In addition to the above variables the DY structure function, denoted 
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by W^, also depends on the heavy flavor masses mi and 777.2. The pointhke cross 
section refers to the reaction 

qi + q2^V ^£1+12, (2.6) 

and is given by (see [2, 3]) 

' N AMz sin^ ew cos^ Ow {Q^ - M|)2 + M|r| ' ^ ' 

For completeness we also give the formula for the 7 — Z interference 

n (n^ Am - 1 1 - 4sin^g^ 1 - M| 

a^zlC^,MzJ-^ 6 sin^^^cos^e^^ (g2_^2)2 + ^2r2- l^-^^J 

In the above expression ^vk denotes the electroweak angle. Furthermore Fz and 
represent the total width of the Z and W boson respectively (sum over all decay 
channels) and = 3 {SU{N), color). The partial widths due to the leptonic decay 
of the Z and W are given by 

_ aMzil + il-Asin'Owf) 
^^-^"^ ~ 48sin2^^cos2^^ ' ^^-^^^ 

In the case of W- and Z-production the total cross section can be obtained using the 
narrow width approximation while integrating (2.5) over Q^, i.e. 

S{Q^-M^). (2.13) 



(g2_M2)2 + M2F2^ MvTv 
The total rates afj (sum over all leptonic and hadronic decay channels) are now given 

by 

^1 - Jf A.-^fl., lmMl/s,Mlm\m% (2.14) 
iV 4 sm Ow cos^ Ow s ■' 

of, - ^^^-ni^wls.Ml^mlml). (2.15) 

iV sm Uw s ■' 

When we consider the reaction where the vector boson decays into a specific lepton 
pair £1, ii aYj in expressions (2.14), (2.15) has to be replaced by 

arr'^'^^arjBiv^eA), (2.16) 



where B(V — > £1^2) stands for the branching ratio 

^ - (2.17) 

Notice that all particles into which the vector bosons decay are taken to be massless 
which is a good approximation since Mz » mf, and the top is too heavy to appear 
in the decay products. 

Further since the electrowcak radiative corrections to sin^ 9w are non negligible it is 
better to replace sin^ 9\y appearing in the denominators of the above expressions by 

where Gp = 1.1667- 10^^ GeV~^ (Fermi constant) whereas in the numerators we have 
put sin^ Ow = 0.2258. 

The above definitions imply that the vector- and axial-vector couplings describing 

the strength of the coupling of the vector bosons to the quarks are hidden in the 

definition for in (2.5). The same also holds for the elements of the Kobayashi- 

Maskawa matrix, denoted by Vq^^q^, which we approximate by retaining the Cabibbo 

angles only and putting the remaining angles and phases equal to zero. 

In the case of the Born process as given by reaction (2.6) the DY-structure function 

becomes 

rV IT/ |2 / / , / ^ 



with Vqj^q^ = 1 in the case V — j, Z. The vector- and axial-vector couplings are equal 
to 



3' 
1 

3' 



g 



.w „.w ^ „w „w ^ 



-r=-7 = ^, a:=a-=--^. (2.20) 

We will now list all order parton cross sections due to heavy flavor production 

contributing to reaction (2.1). 

We start with the two-loop corrections to the Born process 

qi + q2^V, (2.21) 

which are presented by the graphs in Figs. 1 and 2. The first one (Fig. 1) contains 
a heavy quark in the triangular loop and it only contributes to Z-production. The 
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W does not contribute because of charge conservation and the same holds for the 
photon due to charge conjugation (Furry 's theorem). The DY structure function is 
equal to [10]- [12] 



1 /a,^2 

2 V TT 



where Cp and 7/ are the SV (N) color factors: Cf 



2N 



(2.22) 

, Tf = ^ and q — u, d, 

s, c , Q — h, t. The functions Gi and G2 are given in Eqs. (2.8) and (2.9) of [11] 
respectively. Notice that in expression (2.22) one has to sum over h and t in order to 
cancel the Adler-Bell-Jackiw axial anomaly. 

The second two-loop correction to (2.21) is given by the vertex correction in Fig. 2. 
It contains the heavy quark loop contribution to the gluon self energy insertion in 
the vertex graph. The DY structure function becomes 

= - r) IV^,,/ + (aj-)') CfT, \ (^)' F{Q\m% (2.23) 

where F{Q'^,m'^) can be found in (Al). The above expression contributes for all 
vector bosons and Vq^q^ — 1 when V — Z. 

Next we have the one-loop corrections to the two-to-two body processes (see Figs. 3 
and 4) 



g + ?(^) 



q{q) + V. 



(2.24) 
(2.25) 



They all contain the triangle heavy quark loop insertion which only contributes to 
Z-production for the same reasons as mentioned above (2.22). For process (2.24) 
(Fig. 3) the DY structure function reads 



TT 



2 r 



{-2 + 2f (Ji(4m7s) - Ji(4mVQ2))} 



4m 



(j2(4mVs) - J2(4mVg^ 



(2.26) 



where Ji and J2 arc given in Eqn. (2.12) of [11]. 
For process (2.25) (Fig. 4) we have the expression 
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S '^Q ^ / 2 V TT 



9{Q^ -4m^)Hi{s,Q\m'^) 



(2.27) 



with Hi and H2 defined in (2.18) and (2.19) of [11] respectively. Like for the DY 
structure function in (2.22) the above expressions (2.26) and (2.27) have to be summed 



6 



over b and t in order to cancel the axial anomaly. 

The next contributions are given by the two to three body reactions (see Figs. 5 and 
6) 

q{h) + q{k2) ^ V{q) + QM+Q2{p2): (2.28) 

g{ki) + g{k2) ^ V{q) + Q,{p,) + Q2{p2). (2.29) 

Since at this moment experiments indicate that rrit > rrih + Mw [9] Ty-production 
has to be treated in a way which differs from the usual procedure applied to Z- and 
7-production. This is due to the instability of the top quark appearing in the internal 
lines of the Feynman graphs in Figs. 5 and 6 corresponding to the reactions 

q + q t + t (2.30) 

q + q i+t (2.31) 
^W- + b, 

g + g t + t (2.32) 

^W+ + b, 

g + g t + t (2.33) 
^W- + b. 

In this case the internal top quark line cannot be described by an ordinary Feynman 
propagator anymore. Therefore it has to be treated as a resonance which implies that 
the top quark propagator has to be replaced by a Breit- Wigner form. However this 
procedure leads to a violation of gauge invariance which is hard to remedy. In the 
case of the total cross section (2.15) one can resort to an approximation which will 
be presented at the end of this section. 

Starting with Z- and 7-production which is described by the Feynman diagrams in 
Figs. 5 and 6 the DY structure function is given by 



S4 S4 



d(j) / ci^sin^ M-. , (2.34) 
/o Jo •' 

where (i, j) = {q, q) or {g, g). The symbol K^j denotes the color factor which is given 
by Kqq = CpTf and Kgg = Tf. In the definition of Wy an average over the initial 
spins and a sum over the final spins is understood. The Kallen function is defined by 
X{x, y, z) = + y^ + z"^ — 2xy — 2xz — 2yz and the kinematical invariants §4, ti, and 
Ui are given by (see (2.28) and (2.29)) 

S4 = (pi+P2)', ii^{k2-qf, ui^{ki-qf. (2.35) 
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In addition to the above invariants the matrix element squared 



M 



V 



depends on 



other kinematical variables which are analogous to the ones defined in Eqn. (4.2) of 

2 



requires an inten- 



[13]. After having performed the traces the computation of 
sive partial fractioning before the angular integration can be carried out. The angular 
integrals can be found in Appendix C of [13]. Notice that in the definition of 
we have removed all the strong and electroweak coupling constants described by gs, 
e, and g as defined in sections 10.6 and 14.5 of [14]. This also includes the typical fac- 
tors which appear in the vertices like —igs^^Ta, —iev'^^/^, — ^V^iq2 7M('^^ + 75^^) 
— j^^o^lfiiVq + Ibdq) (see Eqn. (2.20)). Here V^^gj denotes the Kobayashi-Maskawa 
matrixelement where only the Cabibbo angle 6c has been put to be unequal to 
zero. All these factors are absorbed in the definition of the pointlike cross sections 
av{Q^, My) as presented in (2.7)-(2.9) except for the Cabibbo angles and the cou- 
plings and (2.20) which remain in W^. In the case the vector boson couples 
to massless quarks they can be factored out too like in (2.19) but for massive quarks 



like the heavy fiavors they remain in 
In the case » (mi 



M 



7712 



DY structure functions W^^ and W^g. 



m) one can obtain analytic expressions for the 
For process (2.28) where the vector boson is 
radiated off the incoming light quark lines (sec Figs. 5a,b) one obtains * 



r, 



Att) 



3 1-f 



1 + 



32 



ln(l 



32 



Inf 



80 



32 



In 



^— ^ y ln2(l - r) - y In T ln(l - r) + InV - — ln(l - r) 



160, . 8^. ,^ 32^, , 448 \ 16 

+— In f - -Li2(l - f) - -C(2) + - (1 - f) 
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41n(l-f)-41nf-^) -^lnf+^(H-f) (^2(1 - f) + ^ In^ f 



(2.36) 

If the parton cross section (2.5) is convoluted by the parton densities in order to 
compute the hadronic cross sections, as will be discussed in the next section, one 
approaches a singularity at f = 1. In this case one cannot neglect the mass m in the 
denominator anymore. The resulting terms which even can go as a power of the type 
In^ Q'^ jrr? will be partially cancelled by similar terms arising in the vertex correction 
(2.23) (see (A3)). They all can be described by an expression proportional to a delta 
function which has to be added to expression (2.36). This expression reads as follows 



V,i2) 



= IK 



9192 I 



C.T, (g)^(l 



*The polylogarithms of the type Li„(a;), Sn,p{x) are defined in [15]. 
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68, Q2 32 
m 



454 



(2.37) 



Furthermore one has to replace in (2.36) the singular terms of the type ln*(l — f)/(l — 
f) by (ln*(l - f)/(l - f))+ with the definition 

When the vector boson is radiated ofi^ the final state (see Figs. 5c,d) one finds in the 
limit » 



(11(12 



2 r 



(l + r) 



--Li2(-f) - -C(2) + - In^ f - - In f ln(l + f) | 



+ ^(3 + 4f + 3f2)lnf+^(l-f2) 



(2.39) 



Finally in the case of Z-production one can also find an asymptotic expression for 
the interference terms between diagrams 5a,b and 5c,d. It is given by [2, 11] 



2 r l + f2 

16- In f + 32f In f + 16(3 - f) 

1 — T 



(2.40) 

For the gluon-gluon fusion process (2.29) (see Fig. 6) one finds the following asymp- 
totic expression 



W^'W (f = IK 



V, 



{8(1 + 4f + 4f2)• 



lnf+16(l-f)(l+3f)}ln'^+| -32(l + 4f + 4f2) f Li2(l-f) + 



Inf ln(l -f) - ^ln=^f ) - 64(1 - f)(l + 3f) ln(l - f) + 



> n2 

8(1 + 8f - 4f2) Inr + 4(1 - f)(7 + 67f) Un + 4(1 + 4f + 4f^) 



16Li3(l - f) - 4Li2(l - f) Inf - 16Li2(l - f) ln(l - f) + 



41n^f ln(l -f) -8W(1 -f) Inf) +4(l + f) (| 8Li2(-f) + 
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Inf ln(l + f ) + 4(1 + ff - 165i,2(-f ) - 16Li2(-f ) ln(l + f ) 



-8C(2)ln(l + f) + 121n2f ln(l + f)-81n2(l + f) Inf - 32(1 + lOf + 



7f2) 5i,2(l - f) - 32(1 + 2f - f')Li3(-f ) - 16(1 + 2f - 2f2)C(3) 
-16(5 + 4f - 14f^)Li2(l - f ) + 32(2 + 4f + f2)Li2(-f) Inf + 16C(2) • 

(3 + lOf + lOf 2) In f + 16(5 + 9f - 12f 2)C(2) - ^(3 + 8f + 8f ^) In^ f 

-8(3 + 7f + 4f 2) In^ f - 64(1 - f )(1 + 3f ) ln^{l - f ) + 16(1 + 8f - 4f 2) 
In f ln(l - f ) - 4(23 + 64f - 105f 2) In f + 8(1 - f ) (7 + 67f ) ln(l - f ) 

-8(1 -f)(16 + 49f) 



Ar2 



- 1 



4(1 + f )2 I 165i,2(-f ) + 24Li3(-f ) + 16C(3) 



16, 



+— Li2(-f ) - 24Li2(-f ) Inf + 16Li2(-f ) ln(l + f ) + 8C(2) ln(l + f ) 

+ -C(2) - 121n2f ln(l + f) + 81n2(l + f) Inf + — Inf ln(l + f) j 
3 3 J 

-32(1 - f )^^i,2(l - r) + ^(-2 + 2f + 25f2) W f - ^(6 + 38f + 75f^) Inf 



-(1 -f)(188 + 764f) 



(2.41) 



As we will see in the next section some of the above approximation turn out to be 
very useful for Z-production accompanied with bb quarks because — M| >> m^. 
As has been discussed below (2.29) 1^-production has to be treated in a different way 
as has been done above for Z- and 7-production. Here one has to make a distinction 
between initial and final state emission of the VT-boson. In the case the W^-boson 
is radiated off from a light quark in the initial state, described by the graphs in 
Figs. 5a,b, the DY structure function is given by W^^ in (2.34). However if the 
VF-boson is the decay product of the top or anti-top quark in the final state like 
in Figs. 5c, d or Fig. 6 one has to resort to different methods. In this paper we are 
only interested in the total cross section. Hence we can follow the same procedure 
as is outlined in (2.13)-(2.15). First we neglect the graphs where the W is emitted 
from the bottom quark because the latter is far off-shell and apply the narrow width 
approximation to the Breit-Wigner form of the top quark in reactions (2.30) and 
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(2.31). This is a reasonable approach because the width of the top Ft — 1.41 GeV^ 
is much smaller than its mass rrit = 174 GeV [9]. Following the above procedure for 
1^-production in quark- antiquark annihilation the total cross section is then given 

by 



<7Y,^<7UQQ^tt)B{t^Wb), 
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with B{t Wh) 1 and cTtot(?? tt) [19, 25] is equal to 

' 2m2^ 



An 1 1 



4m2 



1 + 



(2.42) 



(2.43) 



In the case of the gluon-gluon fusion process (2.29) we proceed in the same way. 
Neglecting the emission of the W from the bottom quark and applying the narrow 
width approximation to reactions (2.32), (2.33) we get 



w 

99 



atot{gg^tf)B{t^Wb), 



where (Ttot{gg ^ tt) is given by [19, 25] 



Am? 



(2.44) 



+ 1 + 



4m^ 8m^ ' 



s P 



1 iv^ r 



'2 lOm^^ 

3 ^ TT 



Arm? 



with 



+ -:^ln|/(s) 



y{s) 



1 + ^l-Am?/s 
1 - Y^l - Am?/s 



(2.45) 



(2.46) 



^Tt is related to mt using the formulas in [18]. 
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3 Hadronic cross sections 



In this section we want to discuss the heavy flavor (6 and t) contribution to the total 
cross section of W- and Z-production at large hadron coUiders. The energies and 
the coUiders under study are given by ^/s — 0.63 TeV {SppS, pp), ^/s — 1.8 TeV 
(TEVATRON, pp) and ^/s = 16TcV (LHC, pp). Investigated will be the part of 
the total cross section constituted by the heavy flavor contribution. We also want to 
know how the latter, which is of order a^, compares with the light parton contribution 
calculated in the same order of perturbation theory which has been studied in the 
past (see [2, 3]). Finally we want to study the validity of the approximation for the 
total cross section oi H1 + H2 —> Z + h-\-h obtained from the formulae in (2.36)-(2.41) 
which are calculated in the limit M| >> m^. 

The hadronic cross section is related to the partonic cross section (2.5) through the 
relation 



I dyi [ dy2 f dz 5{t - yiy-2z) yiy-2 f"^{yu n'^) ff^{y-2, 1^^) 
77 .^0 Jo JO 



dQ- 



:{Q^/yiy2S,Q^,ml,ml,iJ,'^), (3.1) 



where f^iy^lJi^) denotes the density of parton i inside the hadron H which depends 
aside from y also on the factorization (renormalization) scale Notice that for con- 
venience we have put the renormalization scale equal to the factorization scale. 
In the case heavy flavor production is treated in lowest order, as we do in this paper, 
da/dQ"^ is independent of the factorization scale (Born approximation). However 
since this approximation is of order it does depend on the renormalization scale. 
In the case of light partons in the initial and final state one has to perform mass fac- 
torization to dcFij/dQ"^ in order to remove the coUinear divergences and this quantity 
has to be replaced by the DY coefficient function which has been calculated up to 
order in [2, 3]. Therefore in addition to the renormalization scale the latter also 
depends on the factorization scale. 

In our calculations we have chosen the MS-scheme for the coupling constant as as 
well as for the DY-coefRcient function calculated up to order in [2] . For the parton 
densities we have chosen the MS- version of the set MRS(H) [26] with = 230 MeV. 
Furthermore we use the two-loop corrected running coupling constant as with the 
QCD scale A mentioned above. Since we only consider top and bottom quark pro- 
duction we have put the number of light flavors equal to four, i.e. n/ = 4. Finally 
we have set the factorization (renormahzation) scale jj? — where — My. For 
the electrowcak parameters we have taken the following values: Mz = 91.196 GeV, 
Mw = 80.24 GcV, Gp = 1.1667- lO^-'^GeV^^ sin^^vi^ = 0.2258 and sin^^c = 0.0484. 
The masses of the heavy flavors are given by m?, = 5 GeV and mj = 174 GeV. 
To calculate the total cross section for Z-production we have integrated expression 
(3.1) over Q'^ and used the narrow width approximation (2.14). In tables 1-3 we have 
listed the various contributions coming from the partonic subprocesses discussed in 
the previous section and compared them with the light parton cross section corrected 
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up to order a'^. The tables reveal that at smaller energies, i.e. ^/s — 0.63 TeV the 

total heavy flavor cross section is dominated by the vertex correction (2.22) given by 
Fig. 1 and the subprocess q + q Z + b + b (2.34) depicted in Fig. 5. As we will show 
later on the importance of the last process is wholly due to initial state radiation 
(Figs. 5a,b). At larger energies like i/i = 1.8 TeV also the process g + g —>■ Z + b + b 
(2.34) (Fig. 6) becomes important. The latter even overwhelms the other reactions 
when — 16 TeV (LHC). This can be traced back to the gluon density which steeply 
rises at very small r = My/s. Finally we observe that the processes with top quarks 
in the final state are completely unimportant which is due to the limited phase space 
available even for energies as large as ^/s — 16 TeV (LHC). 

The reason that the virtual correction in Fig. 1 plays an important role can be in- 
ferred from the expression in (2.22). Here the first term Gi only contributes in the 
case of the bottom loop (M| > 4ml) whereas the second term G2 only contributes 
for the top loop (M| < 4m|). In [11] it has been shown that for » rn? the 
function Gi vanishes like 

whereas the function G2 has the following asymptotic behavior for » 

which means that this correction is dominated by the top-loop contribution. 
In tables 1-3 we have also listed the results coming from the light parton contribution 
calculated in [2]. We observe that the heavy flavor part of the DY cross section is 
very small and it amounts to 0.2% (v^ = 0.63 TeV), 0.3% (v^ = 1.8 TeV), 1% {^s = 
16 TeV) of the light quark contribution. Even if we compare the heavy flavor part, 
which is an order correction, with the corresponding light parton contribution in 
the same order of perturbation theory one discovers that it is very small except when 
the energy gets very large. It amounts to 5% {^/s = 0.63 TeV), 11% {^/s — 1.8 TeV), 
and 770% (absolute value) {^/s = 16 TeV) of the «^ correction to the cross section 
which is due to light partons. This means that only at LHC energies heavy flavor 
production is more important than the light parton subprocesses contributing in order 

For the PF-cross section we proceed in the same way as for Z-production except that 
here we also have to make the narrow width approximation for the top quark. This 
is needed when the W is radiated off the top-quark in the flnal state exhibited by 
the graphs in Figs. 5c, d and 6. Furthermore the graphs containing the triangle heavy 
flavor loop hke Figs. 1, 3, and 4 do not contribute. The results are given in tables 4-6. 
Like in the case of Z-production the process qi + q2 ^ W + b + b (2.34) (Figs. 5a,b) 
is dominant at lower energies although also the vertex correction (2.23) contributes 
a little bit. When the energy gets larger like in the case of LHC also the process 
g + g ^ t + t with t W'^b and t W~b (2.44) (Fig. 6) becomes important due to 
the gluon density which becomes very large at small r = M^/s. As in Z-production 
the heavy flavors give a small contribution to the VF-cross section. The latter is of the 
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(3.2) 



(3.3) 



same magnitude as in the Z-cross section and it amounts to 0.1% {y/s — 0.63 TeV), 
0.1% {^/s = 1.8 TeV), and 0.2% {^/s = 16TcV) of the hght parton contribution. If 
we make a comparison with the order part of the hght parton contribution these 
numbers become 2%, 6%, and 90% respectively which are however smaUer than in 
the case of Z-production. 

Besides vector boson production we have also studied the heavy flavor part of the 
DY cross section da/dQ^ at small < 60 GeV) where the virtual photon 

dominates the cross section. In Fig. 7 we have plotted the ratio R 

da , , , da 



{u,d,s,c,g) + -—{t,b) 

R{Q') = ^^-^ , (3.4) 

-Xu,d,s,c,g) 



dQ'' 

at three different energies, i.e. ^/s = 0.63, 1.8, 16 TeV. Like in the case of Z- and 
ly-boson production the contribution of the heavy flavors to the DY cross section 
is very small. When the energy increases it grows from 0.1% to 0.5% of the part 
constituted by the light parton contributions. 

Using the MRS(H) [26] parton densities we also calculate the cross sections for Z- and 
W-production where the lepton pair into which the vector boson decays is measured. 
The results are presented in table 7 {SppS, \fS = 0.63 TeV) and table 8 (TEVA- 
TRON, = 1.8 TeV). They are obtained by multiplying the total cross sections 
a^ and a^ in tables 1 to 6 by the branching ratios B{Z — > e+e") = 3.35 • 10~^ and 
B(W ez/g) = 0.109 respectively. Furthermore we have also listed the experimental 
data obtained from the groups UAl [27], UA2 [28] {SppS) and CDF [29] (TEVA- 
TRON). As is expected from the previous discussion the contribution from the heavy 
flavors h and t to the cross section is extremely small when compared with the light 
parton part. 

Before concluding this section we have also studied the cross sections obtained from 
(2.36)-(2.41) which are derived in the limit » m^. In practice these formulae are 
only applicable to the reaction H1 + H2 Z + b + b where M| >> ml. This inequal- 
ity does not apply to top production also when the Z is replaced by the W. In the 
case a photon appears in the intermediate state the above formulae are also not very 
useful because in the region >> the production rate is too low. In table 9 
we have compared the results obtained from the exact cross section represented by 
(2.23), (2.34) with those predicted by the asymptotic formulae in (2.36)-(2.41). First 
the table shows that the whole contribution to the cross section of the subproccss 
g + g— > Z + b + b is given by the initial state radiation of the Z-boson as depicted in 
Figs. 5a,b which also includes the virtual contribution in Fig. 2. This is due to the 
large logarithmic terms In'^ (Q"^ / m'^) which show up in (2.36) and (2.37). Furthermore 
the approximations to the initial state radiation process work rather well in partic- 
ularly when the energy increases. The difference between the exact (2.23), (2.34) 
and the approximate cross sections (2.36), (2.37) amounts to 30% {^/s — 0.63 TeV), 
17% {y/s = 1.8 TeV) and 16% {y/s = 16 TeV) of the exact value. The main reason 
for the difference is that the cubic logarithmic term ln^(Q^/m^), arising in the exact 
expressions (2.23), (2.34) in the limit » m^, will only cancel when is taken 
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to be really small with respect to . This is apparently not the case for = M| 
and w? — ml. Therefore the sum of the two exact expressions is not quite equal to 
the sum of the approximations (2.36) and (2.37) which behaves asymptotically like 
ln^{Q'^ /m^). In practice the asymptotic limit is only reached for < l(GeV/c)^ 
which is an order of magnitude less than = (5GeV/c^)^. A similar problem 
shows up in the contribution originating from the interference between initial state 
(Figs. 5a,b) and final state (Figs. 5c, d) radiation of the vector boson. In spite of 
the fact that here no large logarithms of the type In'' (Q"^ / rrP) appear in the final 
expression (2.40) the various integrals contributing to the interference term contain 
these type of logarithms so that in this case one observes an incomplete cancellation 
too. From table 9 we infer that the approximation to the interference term is not 
so bad when ^/s — 0.63 TeV but it becomes worse at higher energies. Apart from 
the incomplete cancellation mentioned above this is also due to the fact that the 
quality of the approximation depends on the values for s = yiy2S in the partonic 
cross section da/dQ"^ (3-1)- When the energy ^/s increases it may happen that the 
product of the parton densities in (3.1) probe the s-region where the approximation 
fails. Fortunately the interference term gives a negligible contribution to the cross 
section ofg + g^ Z + h + h. The latter also holds for the final state radiation of the 
Z-boson as depicted by the graphs in Figs. 5c, d. However here the approximation 
to the the cross section (2.34) as given by (2.39) becomes excellent which is inde- 
pendent of the energies under consideration. The most simple explanation for this 
phenomenon is that here large logarithms of the type Itl' {Q"^ j m"^) neither appear in 
the final result (2.39) nor in the separate integrals contributing to this expression. 
Therefore in this case no cancellation of large logarithmic terms has to occur. Finally 
table 9 reveals that the approximation to the cross section of the gluon-gluon fusion 
process g -\- g ^ Z -\-h -\-h (Fig. 6) turns out to be very good for all energies under 
consideration. Although the approximate cross section (2.41) behaves quadratically 
in hi{Q'^ /m?) it is closer to the exact result as discovered for the initial state radiation 
of the Z-boson ing + g— > Z + h + h. This is mainly due to the fact that in the former 
process no cancellation of leading logarithmic terms occur. 

In general we can conclude that irrespective of the energies considered the approxi- 
mations work rather well for Z-production accompanied by bottom quarks although 
this statement is more valid for g + g ^ Z + h + h (Fig. 6) than for q + q ^ Z + h + h 
(Fig. 5). Since the large logarithms ln'^((5^/m^) dominate the cross sections for 
\fs > 1.8 TeV the bottom can be treated as a light quark. Further we have also 
studied the limit » w? for the charmed quark cross sections. In this case the 
approximations (2.36)-(2.41) are even better than for bottom production. This even 
holds for i/s = 0.63 TeV. The logarithms of the type ln(Q^/m^) dominate the charm 
and bottom cross sections and they give rise to large corrections. Therefore they 
have to be removed by mass factorization and subsequently to be absorbed, after 
resummation via the renormalization group equations, into the charm and bottom 
densities in the hadron. 

Summarizing the content of this paper we have computed all order a1 contributions 
to the DY cross section which can be attributed to heavy fiavors. All virtual as well 
as radiative processes have been considered. In this way we have extended the work 
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done in [10, 11] where only the contributions characteristic for Z-production have 
been considered. 

From the results obtained in this work one can conclude that the contributions of the 
heavy flavors b and t to the DY cross section in particularly to vector boson produc- 
tion are very small. They are on the one percent level in the case of Z-production 
provided the energy is very large which will only happen when the LHC is put into 
operation. Furthermore we have shown that for ^/s > 1.8 TeV the cross sections 
(2.36)-(2.41) derived in the limit » can be applied to Z-66-production. This 
means that the bottom quark can be treated as a light flavor at large collider ener- 
gies. The heavy flavor cross sections will only become observable if the vector boson 
as well as the heavy quarks are detected. This will happen for the LHC where the 
process > Z + b + b constitutes an important background for Higgs production 

[30]. 
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Appendix A 

In this appendix we present the two-loop vertex correction defined by F{Q'^,m'^) in 
(2.23). It contains the gluon self energy contribution with the heavy flavors appearing 
in the subloop (Fig. 2). The vertex correction reads 

2Li2 - 2C 2 + - ln=^ + 16 — - - 



952 



2 



3355 , 

(Al) 



where 



- = f+4^. (A2) 

Further the following asymptotic expansions turn out to be useful. In the limit 
w? « one obtains the expansion 

mn2 2^ 4 38 Q2 530 \ Q2 



m-*0 



3355 152^, , 16^,^, , , ^, 



In the case that m » Q the vertex correction behaves like 



m"' V 225 45 



m— >cx3 



which shows that heavy flavors decouple from F{Q'^,m?) when the quark mass gets 
infinite. 
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Table 1 

Contributions to the total cross section for Z-production at -^/s — 0.63 TeV 
{as{Mz) = 0.107). 



subprocess 


equation 


(nb) 


q + q^ Z 


(2.22) 


3.13-10-^ 


q + q^ Z 


(2.23) 


1.40 • 10-^ 


q + q^ Z + g 


(2.26) 


1.13-10-^ 


9 + ^ Z + q(q) 


(2.27) 


-1.77-10-^ 


q + q^ Z + b + b 


(2.34) 


1.47-10-=^ 


q + q^ Z + t + t 


(2.34) 


1.35 • 10-^^ 


g + g ^ Z + b + b 


(2.34) 


6.39 • 10-5 


g + g ^ Z + t + t 


(2.34) 


1.25 • 10-16 


a^{b,t) 




4.90 • 10-=^ 


a^{u, d, s, c, g) = 1.54 (Born) + 0.41 (0(q;,)) + 0.10 (0(q;J) = 2.05 



Table 2 

Contributions to the total cross section for Z-production at = 1.8 TeV {as{Mz) — 
0.107). 



subprocess 


equation 


(nb) 


q + q^Z 


(2.22) 


5.33 ■ 10-3 


q + q^Z 


(2.23) 


4.86 ■ 10"^ 


q + q-^ Z + g 


(2.26) 


3.67-10-^ 


g + q{q) ^ Z + q{q) 


(2.27) 


-1.61 ■ 10"^ 


q + q ^ Z + b + b 


(2.34) 


8.57-10-3 


q + q^Z + t + t 


(2.34) 


3.58 ■ 10-6 


g + g Z + b + b 


(2.34) 


4.71 ■ 10-3 


g + g ^ Z + t + i 


(2.34) 


5.50 ■ 10-^ 


a^{b,t) 




1.93-10-2 


a^{u, d, s, c, g) = 5.34 (Born) + 1.05 (0(as)) + 0.17 {0{aj)) = 6.56 
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Table 3 

Contributions to the total cross section for Z-production at a/s — 16TeV {as{Mz) — 
0.107). 



subprocess 


equation 


(7^ (nb) 


q + q^ Z 


(2.22) 


1.21 • 10-2 


q + q^ Z 


(2.23) 


5.02 • 10-3 


q + q^ Z + g 


(2.26) 


8.06 • 10-^ 


9 + (1{q) ^ Z + q{q) 


(2.27) 


-1.57-10-3 


q + q^ Z + b + b 


(2.34) 


9.28 • 10-2 


q + q^Z + f + f 


(2.31) 


3.17 • 10- ^ 


g + g Z + h + h 


(2.31) 


5.85 • 10"^- 


g + g ^ Z + t + i 


(2.34) 


9.23 • 10"^ 


a^ib,t) 




6.95 • 10-1 


a^{u, d, s, c, g) = 55.2 (Born) + 7.45 (0(q;,)) - 0.09 (0(q;^)) = 62.6 



Table 4 

Contributions to the total cross section for + H^" -production at = 0.63 TeV 
{as{Mw) = 0.109). 



sul)])r()cess 


(xiuatioii 


a"' (nl)) 


qi + 1-1 ir 


(2.23) 


1.32 • K)-"' 


gi + ^2 W + h + h 


(2.31) 


1.11 • IC)-'' 


qi + q2^W + t + t 


(2.34) 


1.94 • 10-11 


q + q t{t)+t{t) 

U W+{W-) + b{b) 


(2.42) 


1.43 • 10"^ 


g + g ^ t{t)+t{t) 

U W+{W-) + b{b) 


(2.44) 


1.52 • 10-9 


a^{b,t) 5.43-10-3 


a^{u, d, s, c, g) = 5.09 (Born) + 1.34 (0(a.)) + 0.33 (^(a^)) = 6.76 
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Table 5 

Contributions to the total cross section for W'^ + VF~-production at ^/s — 1.8 TeV 
{as{Mw) = 0.109). 



subprocess 


equation 


(nb) 




(2.23) 


4.67 • 10-=^ 


qi + q2^W + b + b 


(2.34) 


2.51 • 10-2 


qi + q2 + t + t 


(2.34) 


6.36 • 10-6 


q + q t{i)+i{t) 

U W+iW-) + b{b) 


(2.42) 


1.10-10-3 


g + g ^ t{t)+t{t) 

U W+{W-) + b{b) 


(2.44) 


1.28 • 10-^ 






3.10- 10-2 


a^{u, d, s, c, g) = 18.0 (Born) + 3.47 (C(aJ) + 0.50 (C(a^)) = 22.0 



Table 6 

Contributions to the total cross section for + PF" -production at ^/s = 16 TeV 
{as{Mw) = 0.109). 



subprocess 


equation 


(nb) 




(2.23) 


4.79 • 10-2 


qi + q2^W + b + b 


(2.34) 


2.76 • 10-^ 


qi + q2^W + t + t 


(2.34) 


6.26 • 10-^ 


q + q t{t)+t{t) 

U W+iW-) + bib) 


(2.42) 


2.00 • 10-2 


g + g ^ t{t)+t{t) 

U W+{W-) + b{b) 


(2.44) 


1.77-10-^ 


a^{b,t) 




5.22 • 10-^ 


a^{u, d, s, c, g) = 185 (Born) + 24.8 


{0{as)) - 


0.6(O(q;2)) = 209 
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Table 7 

Bz(T^ and Bw for the SppS [27, 28] with Bz = B{Z e+e') = 3.35 • IQ-^ and 
Bw = B{W eue) = 0.109, ^ = 0.63 TeV. 





Bza^ (pb) 


Bw'j'^ (pb) 


UAl 


58.6 ±7.8 ±8.4 


609 ± 41 ± 94 


UA2 


65.6 ±4.0 ±3.8 


682 ± 12 ± 40 


By cr^{u, d, s, c, g) 


68.7 


737 


Bva^{b,t) 


0.164 


0.592 



Table 8 

Bz(T^ and Bw for the TEVATRON [29] with Bz = B{Z e+e') = 3.35 • 10"^ 
and Bw = B{W eVe) = 0.109, ^/s = 1.8 TeV. 





Bza^ (nb) 


Bw (J^ (nb) 


CDF 


0.214 ±0.011 ±0.020 


2.20 ±0.04 ±0.20 


By s, c, g) 


0.220 


2.40 


By a^(b, t) 


6.47 • 10-^ 


3.38 • 10-=^ 
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Table 9 

Comparison of the exact versus approximate cross sections for Z- and 66-production 
(a,(Mz) = 0.107). 









V 


^ = 0.63TeV 












exact 


(nb) 




"app. 


(nb) 






^Z + b + b 


1.58- 


10-3 


(2.23), (2.34) 


1.11 ■ 


10-3 


(2.36), (2.37) 






2.41 ■ 


10-6 


(2.34) 


2.48 ■ 


10-6 


(2.39) 






-9.09 


■ 10-6 


(2.34) 


-1.26 


■ 10-5 


(2.40) 




^Z + b + b 


6.39- 


10-5 


(2.34) 


5.80- 


10-5 


(2.41) 








l.STeV 












exact 


(nb) 




'^app. 


(nb) 




q + q 


^Z + b + b 


8.88- 


10-3 


(2.23), (2.34) 


7.40- 


10-3 


(2.36), (2.37) 






4.31 ■ 


10-5 


(2.34) 


4.35- 


10-5 


(2.39) 






-1.40 


■ 10-5 


(2.34) 


-2.14 


■ 10-5 


(2.40) 


9 + 9 


^Z + b + b 


4.71 • 


10-3 


(2.34) 


4.56- 


10-3 


(2.41) 










/i = 16 TeV 












^exact 


(nb) 






(nb) 




q + q 


^Z + b + b 


9.54- 


10-2 


(2.23), (2.34) 


8.01 ■ 


10-2 


(2.36), (2.37) 






5.93- 


10-^ 


(2.34) 


5.97- 


10-^ 


(2.39) 






-6.88 


• 10-6 


(2.34) 


-2.30 


■ 10-5 


(2.40) 


9 + 9 


^Z + b + b 


5.85- 


10-1 


(2.34) 


5.77- 


10-1 


(2.41) 
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Figure captions 

Fig. 1 Two-loop graph, containing tlie lieavy flavors in the triangular subloop, which 
contributes to the subprocess q + q ^ Z (2.22). 

Fig. 2 Two-loop graphs, containing the heavy flavors in the gluon self energy, which 
contributes to the subprocess q + q ^ V (2.23). 

Fig. 3 One-loop graph with the heavy flavors in the triangle contributing to the 
subprocess q + q ^ Z + g (2.26). 

Fig. 4 One-loop graph with the heavy flavors in the triangle contributing to the 
subprocess g + q{q) Z + q{q) (2.27). 

Fig. 5 Diagrams contributing to the subprocess q + q^V + Qi + Q2; a, b: initial 
state radiation; c, d: flnal state radiation (2.34). 

Fig. 6 Diagrams contributing to the subprocess g + g^V + Qi + Q2 (2.34). 

Fig. 7 The heavy flavor content, represented by R (3.4), of the cross section da/dQ"^ 
for the processes: p+p — > 7* + "X" ^/s = 0.63 TcV (solid hne), ^/s = 1.8 TeV (dotted 
line) and p + p ^ 7* -I- "X" ^/s = IQ TeV (dashed line). 
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